Geometric phase contribution to quantum non-equilibrium many-body dynamics 
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We study the influence of geometry of quantum systems underlying space of states on its quantum 
many-body dynamics. We observe an interplay between dynamical and topological ingredients of 
quantum non-equilibrium dynamics revealed by the geometrical structure of the quantum space of 
states. As a primary example we use the anisotropic XY ring in a transverse magnetic field with an 
additional time-dependent fiux. In particular, if the fiux insertion is slow, non-adiabatic transitions 
in the dynamics are dominated by the dynamical phase. In the opposite limit geometric phase 
strongly affects transition probabilities. We show that this interplay can lead to a non-equilibrium 
phase transition between these two regimes. We also analyze the effect of geometric phase on defect 
generation during crossing a quantum critical point. 



PACS numbers: 

Introduction. — The profound interplay and interrela- 
tion of geometry and physics was the focus in both fields 
since creation of General Theory of Relativity, in which 
quantities responsible for the geometry of space-time are 
determined by the physical properties of the matter living 
in this space and vice versa. The relevant geometry lan- 
guage in this case is a Riemannian geometry. Gauge prin- 
ciple of classical gauge theories found its natural descrip- 
tion and nice interpretation in terms of theory of fiber 
bundles, a subject of differential geometry Jj. Monopoles 
and instantons of the gauge theory have profound topo- 
logical meanings which is the property of defining fiber 
bundle. Many of these notions appeared in various con- 
densed matter systems at equilibrium. Thus, defects in 
He and liquid crystals are classified according to the ho- 
motopy theory, certain phase transitions are associated 
to proliferation of topological defects. Topology plays a 
vital role in e.g. Hall effects and topological insulators. 

Another intriguing phenomena, emerging in quantum 
mechanics, that relates geometry and physics is the Berry 
phase. When a Hamiltonian is adiabatically driven, its 
eigenstates acquire not only the familiar dynamical phase 
factor, but additionally a phase factor that depends only 
on the geometry of the phase space of the Hamiltonian, 
namely the Berry phase It can be observed in inter- 
ference experiments and in the Aharonov-Bohm effect. 
The deep geometrical significance of the Berry phase was 
revealed as well Q. Therefore, it is also referred to as 
the geometric phase. We point that while the Berry 
phase is usually associated with adiabatic processes, the 
geometric phases describe transformations of arbitrary 
eigenstates and are thus not tied to the adiabaticity. In 
condensed matter probably its most transparent man- 
ifestation is in the Haldane phenomena (a presence or 
absence of the excitation gap in ID spin chain depending 
on the value of spins). 

Until now all these manifestations of the geometric 
phase were associated to equilibrium and adiabatic phe- 
nomena. Here we demonstrate for the first time a direct 



relevance of topology and geometry of the quantum space 
of the many body system for the measurable quantities 
defining a non-equilibrium evolution of the system far 
from the adiabatic limit. We show that these effects are 
very significant in the regions close to the quantum phase 
transition. We thus demonstrate a profound interplay of 
geometry and topology of the phase space of the quantum 
many-body system in its out of equilibrium dynamics. 

In equilibrium, a Riemannian structure is introduced 
to quantum mechanics by the Quantum Geometric Ten- 
sor (QGT) The QGT Q^^ is defined for an arbi- 
trary eigenstate \n) by 

Qtiu{\ \n)) := {n\(f^di,\n) - {n\^^\n) {n\dt,\n) , (1) 

for /i, J/ = 1, . . . ,p, labeling the system's parameters 
which form a manifold M. Its real part is a Rieman- 
nian metric tensor g^i, on M that is related to the fi- 
delity susceptibility which describes the systems response 
to a perturbation and therefore is an important quan- 
tity, e.g. in the study of Quantum Phase Transitions 
(QPT) [5[. The imaginary part is related to the 2-form 
(Berry curvature) F^^ := d^Ai, — d^A^ = 23(5^j/, where 
v4^(A, \n)) := j(n|i9^|n) is the connection 1-form. Geo- 
metric phase of the state \n) is given by its integral 
along a closed loop C in parameter space 7„ — A^dA^. 
It is easy to check that after a simple gauge transforma- 
tion, the Schrodinger equation = HItJj) written in the 
instantaneous basis \n) such that 1-0) = J^n '^n\n) can be 
put into the following form: 

On = - ^ Mnm exp [iEnm{t) - «r„„i(t)] a„ , (2) 

where Mnm = {n\dt\m). This equation highlights the 
competition between the dynamical phase Enm(t) = 
Jq [en(T) — e„i(r)] dr and the geometric phase r„m(t) = 
j;;[AA\n))-AA\m))]dr. 

The main purpose of the present work is to demon- 
strate how geometric effects shows up in quantum dy- 
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namics. We do it using an example of a driven XY- 
model which we introduce in the next paragraph. Gen- 
erahzations of some of our results to more generic setups 
are discussed in the Supplementary Information. The 
main findings of our paper is that geometric phase ef- 
fects on transition probabilities are small for slow nearly 
adiabatic driving protocols, i.e. that the leading non- 
adiabatic transitions are determined by the dynami- 
cal phase. Contrary, in the fast limit geometric phase 
strongly affects transitions between different levels. We 
also found that the interplay of geometric and dynami- 
cal phases can lead to non-equilibrium phase transitions 
causing sharp singularities in density of excited quasi- 
particles and pumped energy as a function of the driv- 
ing velocity. This quantum-critical behavior can happen 
without undergoing by the system an actual quantum 
phase transition in the instantaneous basis Q. In the 
limit of slowly driving the system through a quantum 
critical point with an additional rotation in the parame- 
ter space we find that the geometric phase modifies the 
scaling of the observables with the driving velocity and 
enhances non-adiabatic effects. 

The rotated XY spin chain. — Let us consider a stan- 
dard, although rich and illustrative example of XY ring 
in a transverse magnetic field. The Hamiltonian of this 
system 0, S] is defined by 
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with periodic boundary conditions, i.e., ct^^i — a". The 
number of spins N is assumed to be even and the spin 1 /2 
on the site I is represented by the usual Pauli matrices 
af, with a £ z}. Further, the anisotropy for the 

nearest neighbor spin-spin interaction along the x and y 
axis is described by the parameter g and h denotes the 
magnetic field along the z axis. 

At g = this Hamiltonian has an additional U{1) sym- 
metry related to spin-rotations in the XF-plain. At finite 
g this symmetry is broken. Clearly there is a continu- 
ous family of ways breaking this symmetry yielding the 
identical spectrum. The corresponding Hamiltonians are 
related by applying a unitary rotation of all the spins 
around the z axis by angle (p: 

.R(0,z)Ho(5,/^)Rt(0,^), 
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(4) 



with the rotation operator R(0, z) = Jlz^i ^^pC"*'! '''f )■ 
This transformation yields non-trivially complex instan- 
taneous eigenstates, which is a necessary condition for 
existence of the nontrivial geometric phase Q . 

The Hamiltonian Q can be diagonalized using the 
Jordan- Wigner and the Fourier transformations: 



(5) 



with Hfc = {h~ cos pk)a^ + gsmpk(sm2(j)a-'^ — cos 2(/>o'^), 
4 = (e-fe,6l), Pk ^ k = ±l,±2,...,±f and 



Cfc are the Fourier transforms of the fermionic op- 
erators resulting from the Jordan- Wigner transforma- 
tion (see Ref. for details). By applying the Bo- 
goliubov transformation to ([S]) we can map it to a 
free fermionic Hamiltonian with the known spectrum 

<^k{g,h) 



{h^ 



COSPfc)2 + g"^ Sin^pfe. 



The set of quantum critical points of this spin chain are 
determined by the vanishing of the energy gap: 2eko ~ 0' 
where fco is defined by minimizing the excitation energy 
dk^k = 0. This condition defines quantum critical regions 
on A4. For the model (|H) the gap vanishes on the line 
(5 = 0,— l</i<l), marking the anisotropic transition 
and on the two planes S M, /i = ±1), identifying the 
Ising transitions [^, . The anisotropic transition line 
belongs to the Lifshitz universality class since it manifests 
the critical exponents j^i = 1/2 and zi — 2. On the other 
hand the Ising transition planes belong to the d — 2 Ising 
universality class with the critical exponents 1^2 = 1 and 
Z2 — 1 [11]. The points where the critical line and the 
critical plane cross are multicritical points. In Fig. [T] we 
depict the equilibrium phase diagram of the rotated XY 
spin chain in the parameter space (g, h, (j)). 
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Figure 1: The phase diagram of the rotated XY spin chain 
in a transverse magnetic field in cylindrical coordinates: The 
two red planes (ft = ±1) indicate the Ising critical plane, (i.e. 
the associated QPT belongs to the d — 2 Ising universality 
class). Whereas the blue line {g = 0) marks the anisotropic 
transition line. The black bold circle and helix describe the 
two driving protocols we use in this paper. 



Dynamics of the rotated XY spin chain. — We explore 
two driving protocols. The first one is driving the spin 
rotation (f){t) with a constant velocity. This corresponds 
to circular paths in parameter space (see Fig.[T]). It is the 
simplest situation in which a non-trivial geometric phase 
emerges. The second driving protocol consists of driving 
the magnetic field h{t) and the spin rotation 4>{t). This 
results in helical paths in parameter space (Fig. [T]) and 
allows us to study the cross over from the well known 
Landau-Zener scenario (no geometric phase) to the ro- 
tating driving regime (non-trivial geometric phase). 

For either of the protocols we assume (/)(t) = ojt in 
the time interval < t < tf, where a; > is the rate 
of change of the spin rotation. Then the Schrodinger 



3 



equation for the coefficients ai,k and a2.k, that appear 
in the expansion of state in the instantaneous basis, 
— ai.fc|gs)fe + a2,/c|es)i;, becomes a system of hnear 
differential equations with constant coefficients that can 
be solved exactly (see the Supplementary Information). 
From this solution we compute the probability for finding 
the system in the excited state 



Pox.fe = |a2,fc(0/)| = g00(|gs)fc)- 



(6) 



Here f]fc(w) V - ^A^,k] + 4g00(|gs)fc), Ae^ 
Ecs.fc — Cgs.fc = 2efc is the energy difference between 
the excited and ground states of the k-th subspace 
and AA^^k ■= ^^(les)^) - A^{\gs)k) designates the 
corresponding difference of the connection 1-forms: 
A^{\gs)k) = »/c(gs|9</>|gs)fc and A^{\es)k) = ik{es\d^\es)k. 
Further, g00(|gs)fc) is the Riemannian metric tensor of 
the k-th ground state, which also defines the fidelity sus- 
ceptibility along the (j) direction: 

g<^*(|gs)fc) = -fc(gs|90|es)fefc(es|a0|gs)fe = \k{es\d^\gs}k\^ ■ 

With this the total density of excited quasi-particles and 
the energy density of excitations of the entire spin chain 
in the thermodynamic limit can be calculated by 



dk 
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Before proceeding with the detailed analysis of these 
two quantities let us make some qualitative remarks on 
Eq. (i) For the quench of infinitesimal amplitude 

0j — )■ both geometric and dynamical phases are not 
important and the transition probability is simply given 
by the product of the square of the quench amplitude 
and the fidelity susceptibility in agreement with gen- 
eral results [12] ■ (ii) In the slow limit oj <C Aek and 
fixed 0/^1 the geometric phase is still not impor- 
tant while the dynamical phase suppresses the transi- 
tions between levels such that Pcx,k oc g00(|gs)fc)cj^/e^,. 
This result is again in perfect agreement with the gen- 
eral prediction for linear quenches in the absence of geo- 
metric phase 12| given that in this case uj is the velocity 
of the quench, (iii) The most interesting and nontriv- 
ial situation where the geometric phase strongly affects 
the dynamics occurs when both the rotation frequency 
and rotation angle are not small: uj > Ae^, > 1. 
In particular, in the limit — > cx) and (jjf — ^'t- 
recover Pex,fc = 0. This trivial physical fact that in- 
finitely fast rotation can not cause transitions between 
levels actually comes from the mathematical identity: 
(AA^,fe)2 +4g^^(|gs)fc) = [Tr(9^)]2 = 4. For large but 
finite UJ and (f)f ~ 7rn, we find 



(,(|gs)fc)sin^ 
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If the rotation angle is not large n ~ 1 we see that the 
transition probability in this case is directly proportional 
to the square of the product of the geometric and dynam- 
ical phase differences between the ground and excited 
states: 



Pox,fc ~ g00(|gS>fe) 



AEkAj^^, 
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where A^^^k = AA^^k4>f = j^* A^^kdcf) and AEk = 
AcfcT; T — 271 /uj is the rotation period. In the limit 
of large rotation angle at fixed frequency AEk ^ 1 and 
4>fAEk ^ 1 the expression for the transition probabil- 
ity saturates at a value independent of the geometric 
and dynamical phases: Pcx ~ g00(|gs)fe)/2. Interestingly 
this probability is entirely determined by the Riemannian 
metric tensor, i.e. has a purely geometric interpretation. 

From the discussion above we see that if we focus on 
the limit of large 0/ and analyze the transition proba- 
bility as a function of uj we expect a smooth crossover 
between two simple regimes both independent of the geo- 
metric phase: pcx,k ~ g,j„p{\gs) k)uj^ / Ael at a; < Aek and 
Pcx,k ~ g^^{\gs)k)/2 at w > Ack- A similar crossover 
between fast and slow regimes is expected in the many- 
particle situation. Thus one can naively expect that the 
influence of the geometric phase on the dynamics in the 
limit of large 0/ is quite limited. The reality turns out to 
be much more interesting though as we illustrate below. 
In this limit we can simplify the k integrals in Eqs. ([7]) 
using the stationary phase approximation. Then we 
find that the resulting behavior of ricx and Eox exhibits 
a "cusp" at a critical driving velocity uJc determined by 
Wc = 1 — h. This is illustrated in Fig. [2] Because this 



h=0. 




Figure 2: A plot of the total density of excited quasi-particles 
(red line) and the energy density of excitations (blue line) as 
a function of the rotation frequency lj in the limit of many 
rotations </)/ S> 1 and for vanishing magnetic field h = 0, with 
different anisotropy g = 0.01,0.1, 1.0, 10, showing a cusp at 
the critical driving velocity cjc = 1, which can be interpreted 
as a "dynamical quantum phase transition" . 
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singularity is recovered via the stationary phase method 
we expect it to be valid for a class of models with similar 
Hamiltonians. This cusp and the associated "dynamical 
quantum phase transition" is directly related to the 
effect of geometric phase. To understand this let us 
apply a unitary transform tJk{t) ~ diag(e'^"^*^*\ e^"^*^*^), 
to go into a rotating frame, where the geometric 
phase is removed from the Hamiltonian. The resulting 
Hamiltonian in the rotating frame reads rot — 
[{h — cospk) + dt(j)] (7^ + g sinpfe (sin 2(j>a^ — cos 20(7^), 
where the spectrum takes the following form 

Efc.rot — \J {h + u) — cospkY + s\T? pk- From the 
spectrum we see that the Hamiltonian in the rotated 
frame has a quantum phase transition at + Wc ~ 1- 
This transition gives raise to the cusp in Fig. [21 We 
note though that the emergence of the cusp is non- 
trivial since by quenching rotation frequency we are 
pumping finite energy density to the system. In the 
equilibrium this model does not have any singularities 
at finite temperature. Thus this singularity is a purely 
non-equilibrium phenomenon. Further, Fig. [5] illustrates 
nicely that for a small g the regime where riox and Cox 
saturate with uj is close to Wc- However for 5 > 1 the 
dependence of the saturation point on g is approximated 
numerically as uJsAq, h = 0) = 5l.7g°-^^ + 2l.8g^-^'^. 

Another possibility to analyze the interplay of geomet- 
ric and dynamical phases on excess energy and density 
of excitations is to consider the following helical driv- 
ing protocol: {h{t) = St,(f){t) = ujSt), beginning in the 
ground state at = and stopping at t/ = |, i.e. cross- 
ing a quantum critical point. Now 5 plays the role of 
driving velocity, both in h and, for cj ^ 0, in circular 
directions and uj determines the hclicity of the path. For 

0.45| : ^ i i i i 1 




Figure 3: Landau-Zener to Helix: The density of excited 
quasi-particles is plotted as a function S the driving velocity 
for different helicities oj = 0, 0.9, 3, 5, 7, 10, 12 (from down to 
the top) with A'' — 300 spins and an anisotropy of g — 0.9. 
For non zero cu we observe linear scaling regime with S. 



w = we realize the usual Landau-Zener protocol and 



for a; > we describe a helical path in the parameter 
space. In Fig. |3]we present the density of excitations ob- 
tained from an exact numerical integration of the time- 
dependent Schrodinger equation. We recover (Fig. ^ for 
a; = (lowest curve) the scaling riex v^, as expected 
by the Kibble- Zurek scaling argument [l^jQ]- With in- 
creasing w the density of excitations makes a crossover 
to a different linear scaling regime with 6. However, in 
accord with our general discussion in the strict adiabatic 
limit we always observe rtcx oc 

Conclusion. — In summary, we addressed how the 
geometric phase infiuences quantum many-body non- 
equilibrium dynamics. We showed that at intermediate 
of fast driving regimes geometric phase strongly affects 
transition probabilities between levels. We showed that a 
dynamical quantum phase transition can emerge as a re- 
sult of a competition between the geometric and dynam- 
ical phases. This transition manifests itself in the "cusp" 
in the driving velocity dependence of various observables 
(like e.g. the density of excitation and the energy den- 
sity) at finite energy. This allows us to probe quantum 
criticalities "from a distance" , without actually crossing 
them. Such a possibility should be attractive from an 
experimental point of view since the system doesn't need 
to undergo a QPT. We also found that the geometric 
phase modifies the scaling with the driving velocity as 
compared to the LZ scaling. This can be related to ef- 
fective topology-induced interaction between the defects. 
This effect is stronger in the gap-less regions of the phase 
diagram. We also note that our results rely only on the 
geometry of the phase space and thus rather generic. We 
expect that they extend to other protocols where one 
applies a time-dependent unitary transformation to the 
Hamiltonian or other transformation which involves non- 
trivial geometric phase. In particular, similar considera- 
tions apply to the Dicke model realized in Ref. 15]. This 
and possible other generalizations of our results (e.g. for 
open [l6| or turbulent systems) will be discussed in 
a separate work. 
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Supplementary Information 
Effective Hamiltonian in the rotating frame 

Let us consider a more general setup, which extends 
the example studied in the main text. Suppose that the 
system is described by some interacting Hamiltonian Hq, 
which is rotationally invariant with respect to some vec- 
tor coupling A This coupling can represent, for ex- 
ample, an external magnetic or electric field, anisotropic 
interaction constant, couple to a nematic order param- 
eter etc. It can also break an internal U(l) symmetry 
of "Ho like the mixing symmetry between different spin 
components. At nonzero value of A the rotational sym- 
metry is thus broken. Now let us consider a dynamical 
process where this coupling uniformly rotates at a fixed 
magnitude. Then the time dependent Hamiltonian reads 



H{t)^u-\t)n{\o)u{t), 



(10) 



where U (t) is the unitary operator corresponding to this 
rotation. In the rotating frame \ilj{t)) = U\%l}{t)) the ef- 
fective Hamiltonian in the Schrodinger equation picks up 
an additional "centrifugal" term: 



ihdt\il}) = 'HcfflV'), 
Hcff = H(Ao) - 



(11) 
(12) 



where ij) is the rotational angle and a; = is the fre- 
quency. The centrifugal term in the Hamiltonian has a 
number of interesting properties, (i) It is proportional 
to the frequency uj and thus can be used to continuously 
modify the effective Hamiltonian. (ii) At constant fre- 
quency the effective Hamiltonian is time independent, 
(iii) The diagonal components of the centrifugal term are 
given by the connection 1-form of the corresponding en- 
ergy levels. Let us point that the time independence 
of the centrifugal term, which follows its rotational in- 
variance, and its locality imply that the rotations in the 
generic parameter space do not lead to the continuous 
heating even in ergodic non-integrable systems. This sit- 
uation is opposite to e.g. Floquet Hamiltonians. Which 



are usually non-local and lead to constant energy absorp- 
tion in generic interacting systems. Thus we expect that 
the qualitative results of the Paper are valid even if we 
add arbitrary interactions to the Hamiltonian, which pre- 
serve the rotational symmetry but break its integrability. 



Derivation of Eq. ^ 

Here we show how Eq. ([SJ in the main text can be ob- 
tained for a generic two-level system using only the min- 
imal assumption that time dependence enters through 
some U(l) rotation with constant frequency. For such 
a system the Schrodinger equation in the instantaneous 
basis {|gs), jes)} can be written as 

dtCLgs = -(gsjStjes) exp [lE^s.cS) ~ *rgs,cs(i)] Oos, (13) 
dthcs = -(esjStjgs) exp [i£'cs,gs(i) - iFcs,gs(i)] Ogs, (14) 

where Enm,{t) and F„m(i) are the dynamical and ge- 
ometric phases as defined in the main text. No- 
tice that we applied the gauge transformation a„ = 



On exp 



to the coefficients in- 



troduced by: IV'(i)) = ags(i)|gs(i)) + acs(i)|es(t)). In 
agreement with the discussion in the previous section 
all matrix elements appearing in the Schrodinger equa- 
tion ([T^ , ([H]) are time independent for the angle linearly 
changing in time (f>{t) — cot. It is convenient to change 
variables from time t to the angle cj) in (IT^ . P^ : 



d, 



js jes) exp 



(esja^jgs) exp 



UJ 



(15) 
(16) 



differentiating Eq. (1161) one more time with respect to (f) 
and eliminating Ogg, d^dgs gives a second order differen- 
tial equation with constant coefficients: 



Ae 



AA^ d^a^s + g00(|gs))acs = 0. (17) 



From this it is trivial to derive Eq.([5]) from the main text: 



,2 ri 



WMf)\^ = g00(|gs))- 



[if](a.)]' 



(18) 



with n{Lo) := ^[^-AA^]'+4g^^(|gs)). 



[1] In general the Hamiltonian T-Lq can be invariant with re- 
spect to an arbitrary continuous group, not necessarily 
rotations 



